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Large-scale quantum computers will require the ability to apply long sequences of entangling gates to many 
qubits. In a photonic architecture, where single-qubit gates can be performed easily and precisely, the 
application of consecutive two-qubit entangling gates has been a significant obstacle. Here, we demonstrate 
a two-qubit photonic quantum processor that implements two consecutive CNOT gates on the same pair of 
polarisation- encoded qubits. To demonstrate the flexibility of our system, we implement various instances 
of the quantum algorithm for solving of systems of linear equations. 

Quantum computers have attracted tremendous interest because they promise to outperform classical 
computers at a number of tasks\ The long-term goal is to build a universal quantum processor, which 
can be programmed to solve arbitrary problems^. Doing so will require decomposing quantum compu- 
tations into sequences of quantum gates acting on only a few qubits \ In particular, any computation can be 
decomposed using single-qubit gates and a particular two-qubit entangling gate such as controlled-NOT 
(CNOT). In linear- optical quantum computing, two-qubit gates are difficult to implement because photons 
do not interact directly, meaning that the difficulty of a computation can be described by the number of 
entangling gates it requires. 

Entangling gates can be implemented optically by introducing an effective nonlinearity via measurements^, but 
this implies, in the simplest case, that output photons of a two-qubit gate need to be measured to ensure correct 
gate operation, making another application of gates to the same pair of qubits impossible^"^\ Here, we dem- 
onstrate two CNOT gates acting on the same pair of qubits. We implement the two CNOTs differently: the success 
of the first gate is heralded by the successful detection of two ancillary photons^^'^^, while the second gate is 
probabilistic^"^. In distinction to previous work^^ our scheme does not involve qubit recycling which, in principle, 
can be performed deterministically. 

We then use our setup to determine the solution of several two-dimensional systems of linear equations. 
Whereas for a classical computer the number of computational steps needed to solve a linear system is at least 
proportional to the number of variables, a recently developed quantum algorithm could, in some circumstances, 
make the computational time proportional to only the logarithm of the number of variables An important 
difference is that the quantum algorithm calculates the expectation value of an operator associated with the 
solution rather than the solution itself. We note that this algorithm has been implemented previously in photonic 
systems^^, and also using NMR techniques^^'^^. We use the algorithm as an example of the correct operation of two 
consecutive CNOTs, a flexible approach that differs from previous implementations. As we note below, our 
method is as powerful as that in^^. Our implementation thus presents a general quantum circuit with two 
consecutive CNOT gates which can be flexibly adapted to solve other problems. This is a step towards a two- 
qubit universal quantum processor, which would require three consecutive CNOT gates. 

Results 

In our experiment, we use polarisation-encoded photonic qubits (Fig. 1, Fig. 2, and 3b), where |0) and 1 1) denote 
horizontal and vertical polarisation, respectively^^'^^. We produce entangled photon pairs using non-collinear 
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Figure 1 | Experimental setup. Shown is the experimental implementation of two concatenated CNOT gates. The input is set by a polariser, which 
can be followed by a local unitary operation (LU). The two gates are connected by fibers. Different matrices A can be implemented by adapting the LUs, 
and different states I b) by adapting the input state. The figure shows the most general case of two concatenated CNOT gates, combined with general LUs. 
For the implementation of the algorithm, we chose some LUs to be the identity and obtained the case shown in Fig. 3a. 



type-II spontaneous parametric down- conversion (SPDC)^°. A 
pulsed UV laser beam (394.5 nm, 76 MHz, t = 200 fs) passes 
through a 2 mm-thick j5-barium borate crystal, gets reflected, and 
passes through the crystal a second time (see Fig. 1). 

The photons created during the first pass of the laser beam enter 
the first CNOT as the input (control and target) qubits. The state of 
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Figure 2 | Example of a result of a computation on our two-qubit 
quantum processor, (a) The input to the circuit is chosen to be the state 
1 1)10). Then we apply a CNOT gate, followed by a Hadamard gate on the 
control qubit. Finally, we implement the second CNOT gate and 
characterise the output state via two-qubit state tomography, (b) the ideal 
output state and (c) the reconstructed density matrix (real parts), which 
has a fidelity of 0.717 ± 0.030 with the ideal state. The imaginary part is 
below 0.1 and hence not shown. 



these input qubits can be modified using additional local gates and 
polarisers, which in principle allows for the creation of arbitrary 
input states. The first CNOT uses an entangled ancillary photon pair 
and a successful measurement of two ancilla modes heralds that the 
gate has worked correctly on the input qubits^^'^^. The photons cre- 
ated during the second pass act as the entangled ancillary photon 

pairs, in the state + ) = (|00) + | ll))y/\/2. The photons interfere at 

the polarising beam splitters (PBS) as shown in Fig. 1. The PBS 
on the control (target) side is aligned to act in the basis 

{|0),|1)} (^||±) = (|0)±|l))y/\/2|^. The photons are filtered 

spectrally and spatially using narrow-band filters {AX = 3 nm) 
and by coupling them into single-mode fibers. A coincidence detec- 
tion of the ancilla qubits in detectors 3 and 4 in the state |-)3|1)4 
signals a successful gate operation^"^. 

The output photons of the first CNOT (modes I and II in Fig. 1) are 
then guided to the second CNOT. This gate is implemented in a 
destructive way, where a coincident measurement of the output 
photons signals the correct gate operation. The basic element of this 
destructive CNOT is a polarisation-dependent beam splitter (PDBS) 
which has a different transmission coefficient T for horizontally and 
vertically polarised light {Th = I, Ty = 1/3)^"^. If two vertically- 
polarised photons are refiected at this PDBS, they acquire a phase 
shift of 71. Subsequent PDBSs with opposite splitting ratios equalise 
the output amplitudes. This setup, in combination with two half- 
wave plates (HWPs) (Fig. 3b) implements a destructive CNOT. The 
success of this second CNOT is determined by postselection on a 
coincidence detection in outputs 1 and 2. 

Combining these photonic CNOT gates with local unitary opera- 
tions allows us to implement a two-qubit quantum processor. In our 
setup, we implement these local unitary operations using quarter- 
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(a) Simplified circuit for the algorithm for solving specific systems of linear equations. Given a Hermitian matrix A and input I b)y outputs 

bit is measured to be 1. The local unitary R diagonalises A, A = i?^ ^ ^ 

where /li,2 are the eigenvalues of A. (b) Experimental implementation of the circuit in a. The local unitary operations are implemented with the help of a 
combination of half- wave plates and quarter-wave plates, (c) The figure shows a system of linear equations, described by the matrix Ai as well as the 
state vector I bi). (d) The reconstructed density matrix of the experimentally obtained output state \x) is shown. This density matrix is obtained by 

^11 \ /3 \ . 13.. 

n; 1 -i^^ f - 71 1 . We choose the eigenvalues to he Xi = - and ^2 = - by implementing RyiO) correspondingly. The 

fidelity of the reconstructed density matrix is 0.953 ± 0.026. The wire frame shows the theoretical prediction, (e) The quantum algorithm is based on 
determining the expectation value {x\M\x) of some operator Mwith respect to the output state Ix). Therefore, we also show the experimentally 
determined (blue) and theoretical (black) expectation values of several operators M. We choose the operator Mto be the projection on the states 1 0), I +), 
and I respectively, with | + ) = (|0) + |1>)/ \/2, and | +/) = (|0) + /|1)) / Vl. Another example is shown in the SI (Fig. S3). 



choosing the local operation Ri = R^ 



and half- wave plates (QWP and HWP), which can be easily adjusted 
to program the quantum computer (Fig. 1 and 3b). Fig. 2 shows a 
sample run of our experiment; details about experimental para- 
meters can be found in Methods and in the Supplementary 
Information (SI). 

We use our two-qubit quantum processor to solve systems of 
linear equations, which, given a matrix A and a vector b, means 
finding the vector x such that Ax = b. If we rescale the vectors to 
ll^ll = \ \x\ \ = 1 and represent them as quantum states \ b) and \x), the 
task becomes finding \x) such that 



A\x) = \b). 



(1) 



The solution is \x) = A'' \b)/\\A-' \b)\\. 

The quantum algorithm^^ consists of three steps, described in 
Methods. There, we also discuss the simplest possible realisation of 
the algorithm, for a 2 X 2 matrix A, and its reduction to the simplified 
circuit shown in Fig. 3a. This circuit depends on the eigenvalues, 
and ).2y of the matrix A, the unitary R that diagonalises it, 

(2) 



A = R 



0 h 



and the input state \b). 

We have implemented various instances of the algorithm, varying 
all the parameters in A, i^, and | b). The control and target qubits are 
prepared in the states R\b) and respectively, absorbing the local 
operation R into the preparation of the input state. The matrix A can 
be modified by tuning the local operation R (in which case the eigen- 
values stay the same) or by adapting the eigenvalues of the matrix, Xi 
and /I2, which is done by adapting the local unitary operations before 
and after the second CNOT gate (see Fig. 3). Finally, the detection of 
the ancilla qubit in the state 1 1) announces a successful run of the 
algorithm and the preparation of the output qubit in the state \x). The 



success probability is (^1/^2)^' Because we are able to distinguish two 
eigenvalues, our method, although different and using logical fewer 
qubits, is as powerful as that of Cai et aP. 

The output state is analysed using HWPs, QWPs and polarising 
beam splitters, and a full state tomography is performed. Error esti- 
mates are obtained from a Monte Carlo routine assuming Poissonian 
counting statistics. These indicate a lower bound for the actual error 
that takes all the experimental imperfections into account. 

Detailed results of a sample run of the algorithm are shown in 
Fig. 3. Additional results are presented in Fig. 4, where the input 
states and the eigenvalues of A are varied. It can be seen there that 
the performance of the algorithm depends on the input state R\ b). To 
analyse this behaviour, we chose two different input states, | ^1) = 1 1) 
and \b^ = |+) for each A, while keeping R equal to the identity 
matrix. We achieve fidelities of up to 0.981 ± 0.009 for \bi), and 
0.832 ± 0.031 for |^2)- The difference in fidelities arises due to the 
influence of higher-order emissions, as discussed in the SI. Data for 
additional choices of A, i^, and | b) are shown in the SI. The reported 
density matrices contain the raw measured data and have not been 
modified in the post processing. 

Discussion 

We present a two-qubit quantum circuit with two concatenated 
entangling gates acting on the same pair of photonic qubits and 
employ our system to implement the simplest case of the quantum 
algorithm for solving systems of linear equations. With dem- 
onstrating the consecutive application of two CNOT gates on 
the same pair of qubits, we have addressed a technological chal- 
lenge in photonic quantum computation. This development is a 
step towards universal multi- qubit quantum processors using 
single photons as a versatile platform for computation and 
simulation. 
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Figure 4 | The figure shows the solution of the system of hnear equation for matrices with different eigenvalues. Experimentally, these are obtained by 
implementing different values of 6. For all matrices we run the algorithm for two input states I ^i) = II) = (0,1) and |^2) = |+ ) = (l5l)/ V2' We 
achieve fidelities of 0.957 ± 0.010, 0.961 ± 0.013,0.981 ± 0.009 for the input state \ bi) (upper row from left to right) and 0.778 ± 0.031,0.773 ± 0.027, 
0.832 ± 0.031 for the input state I ^2) (lower row, left to right). 



We anticipate that increasing technological capabilities, includ- 
ing the implementation of more than two consecutive CNOTs, 
will allow the extension of the algorithm both to larger systems 
and to more precise phase estimation. Eventually, this will allow 
not only solving larger systems of hnear equations, but also 
equally important algorithms that use it as a subroutine, including 
quantum algorithms for solving nonlinear differential equations^^ 
and quantum data fitting^^. 

Methods 

Theory. We assume, without loss of generality^^, that A is an AT X AT 
Hermitian matrix with eigenbasis {\uj)} and eigenvalues {/Ij}, rescaled so that 



0 < Ij < 1. If I ^> is expanded as \b) = prepare, up to 

normalization, 

j=l 

The quantum algorithm^^ consists of three steps. The most general circuit 
involves phase estimation, a controlled Ry rotation, and a reverse phase 
estimation. Our work differs from the original proposaP"^ in several modifications 
that are needed to implement the algorithm with a limited number of qubits. 
Here, we briefly outline the algorithm while noting our modifications. 

The first step, phase estimation, is a general procedure for decomposing quantum 
states into a particular basis^'^^. It requires an additional "eigenvalue" register of m 
qubits, each initialised to |0>. Phase estimation transforms |^>|0>®™ into 
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Figure 5 | The simplest case of the quantum algorithm for solving systems of linear equations, (a) The complete circuit, as described in Methods, 
with U = exp {2ni 2"~^ A) and 6 = —2 arccos(/li//l2), where ^1^2 are the eigenvalues of A and the integer n depends on the eigenvalues, (b) For 
the algorithm to work perfectly with one eigenvalue qubit, 2 1,2 must be distinguishable by reference to a single binary digit. In that case, 

^ ^ R^ZRy a simplification reflected in this circuit. The middle qubit can now be completely removed. Although controlled 

we follow a different approach, decomposing the controlled rotation to give the 



R-- 



^0 -1 

single- qubit rotations have been implemented in linear optics' 
final circuit shown in Fig. 3a. 



SCIENTIFIC REPORTS | 4:6115 | DOI: 1 0.1 038/srep061 1 5 



4 



y^^^ ^ f^j I l-^;)' where the eigenvalues are stored in the eigenvalue register to a 
precision of m bits (binary digits). 

We restrict ourselves to the simplest case, which involves one state qubit and one 
eigenvalue qubit. The single state qubit implies that |fe> is a two-vector and A a 2 X 2 
matrix, while the single eigenvalue qubit means that only a single bit of the eigen- 
values is computed by the phase estimation. This means that for the algorithm to work 
perfectly, it must be possible to distinguish the two eigenvalues with reference to a 
single bit. Consequently, we choose the two eigenvalues to be of the form O.aO and 
0. a 1, where a is a sequence of bits. The position of the differing bit is what governs the 
choice of the phase estimation unitary U. In particular, if we wish to read out the nth 
bit, we must choose U = exp {2ni A). We note that choosing eigenvalues that 
could not be distinguished using a single bit would lead to additional rounding errors, 
which would, in principle, be avoidable with a larger quantum computer. 

The two-qubit phase estimation transforms the state and eigenvalue qubits, initi- 
alised to I ^>s| 0>£ to y^ ^. _ ^ I ) I ) , where is the n\h bit of Xy Each eigenstate in 
the state register is now entangled to the nSS\ bit of its eigenvalue. 

The second step is the nonunitary map „ ) ^ ) • For this, we introduce an 
additional "ancilla" qubit initialised to |1>, and controlled on the value |%j„) in the 
eigenvalue register, we implement a controlled-i^^(^?j) = exp{ — i6jay/2) rotation on 
the ancilla qubit, where ay is the Pauli matrix. With 9j = — 2arccos (^C2,j~^^ (where C 
< min^ the overall state becomes 



E/*il»i)l^i.«>(^i-c^V 

;=1 



) + ari|i>). 



(4) 



The third step is to run phase estimation in reverse to uncompute giving 

2 

j]/!>,)|0)(^l-C21-^|0) + a-'|l)). (5) 

;=1 

We measure the ancilla, and if we observe a 1 we will have prepared \x) in the state 
register. If we know the eigenvalues, we can maximise the success probability by 
choosing the largest possible C, C = min^- 

The procedure described above involves three qubits, one each for the state, 
eigenvalue, and ancilla registers. The corresponding circuit is shown in Fig. 5, which 
also describes how it can be simplified to the final circuit shown in Fig. 3 a. 

Experiment. For the creation of entangled photon pairs, a mode-locked Coherent 
Mira HP Ti:Sa oscillator is pumped by a Coherent Inc. Verdi V-10 laser. The pulsed- 
laser output (t = 200 fs, /I = 789 nm, 76 MHz) is frequency- doubled using a 2 mm- 
thick lithium triborate (LBO) crystal, resulting in UV pulses of 0.75 W cw average. We 
achieve a stable source of UV pulses by translating the LBO to avoid optical damage to 
the anti- reflection coating of the crystal. Our detectors have quantum efficiencies of 
0.4 and dark count rates between 500 and 3500 counts per second. In our experiment, 
typical visibilities (fidelities) of the emitted Bell pairs of about 0.9 (0.8) and higher- 
order emissions degrade the quality of our gate operations. The back- reflecting 
mirror is continuously moved back-and-forth to avoid any phase correlations 
between of the signal and the noise originating from higher-order photon emissions. 
Additionally, imperfect visibilities on the order of 0.85 to 0.9 of the quantum 
interference (caused by partial distinguishability of the photons) at the PBSs in the 
first gate and the PDBS in the second gate contribute to errors. The measured 
production rate of two-fold coincidences is 18000 s~^ in the forward direction and 
10000 s" ^ in the backward direction (measured after the first gate) . The total fourfold- 
coincidence rate is 1700 h"^ (55 h"^) after the first gate (both gates). 
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